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ABSTRACT 

The impact of strong magnetic fields B > 10 13 G on the thermal evolution of neutron stars is 
investigated, including crustal heating by magnetic field decay. For this purpose, we perform 2D 
cooling simulations with anisotropic thermal conductivity considering all relevant neutrino emission 
processes for realistic neutron stars. The standard cooling models of neutron stars are called into 
question by showing that the magnetic field has relevant (and in many cases dominant) effects on 
the thermal evolution. The presence of the magnetic field significantly affects the thermal surface 
distribution and the cooling history of these objects during both, the early neutrino cooling era and 
the late photon cooling era. The minimal cooling scenario is thus more complex than generally 
assumed. A consistent magneto-thermal evolution of magnetized neutron stars is needed to explain 
the observations. 

Subject headings: Stars: neutron - Stars: magnetic fields - Radiation mechanisms: thermal 



It has been long hoped that the comparison of theoreti- 
cal models for the cooling of neutron stars (NSs) with the 
direct observation of their thermal emission would help to 
unveil the physica l conditions in the interior of these fas- 
cinat ing objects ([Page et al.l [2004: Yak ovlev fc Pethickl 
2004) . Our knowledge of the cooling history of a NS has 
been improving as we were refining the physical ingre- 
dients that play a key role on the thermal evolution of 
NSs. In the past, the field has become more exciting ev- 
ery time that a new relevant idea was introduced (direct 
Urea, superfluidity in dense matter, fast processes due 
to exotic matter, etc.). However, despite the fact that a 
number of NSs are known to have large magnetic fields, 
most studies assumed weak magnetic fields. The main 
reason for this simplification was that the observed dis- 
tribution of magnetic fields in radio-pulsars peaks in a 
region where its effect were thought not to be relevant. 

The increasing evidence that most of the nearby NSs 
with reported thermal emission in the x-ray band of the 
electromagnetic spectru m have anisotropic surface tem- 
perature distributions (|Zavlinl [20071 lHaberll 120071 ) . the 
striking appearance of magnetars (jKaspil 120071 ). and the 
discover y of thermal emission from some high field radio- 
pulsars (|Gonzalez et al.l [20051 ). are indicating that most 
NSs which can be potentially used to contrast theoret- 
ical cooling curves have actually large magnetic fields 
(B > 10 13 G). The conclusion is that a realistic NS cool- 
ing model must not avoid the inclusion of high magnetic 
fields. 

The so-called minimal cooling scena rio (see e.g. 
iPaee et alJl200l lYakovlev fc Pet hick 20Ql for recent re- 
views) defines the cooling model in which the emissivity 
is given by slow processes in the core, such as modi- 
fied Urea and nuclcon-nuclcon Brcmsstrahlung, and en- 
hanced by the neutrino emission from the formation and 
breaking of Cooper-pairs of superfluid neutrons. On the 
other hand, if fast neutrino processes (i.e. direct URCA) 
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take place, the evolution of a NS changes dramatically, 
resulting in the enhanced or fast cooling scenario. Never- 
theless, direct URCA only operates in the inner core of 
high mass NSs for some equations of state. 

In this letter, we want to revisit the minimal cooling 
model considering the effects of magnetic field. If a min- 
imal model must include the minimum number of ingre- 
dients (but all the necessary ones) to explain the obser- 
vations, the magnetic field should be taken into account 
as well. 

The effect of the magnetic field on the surface temper- 
ature distribution caused by the anisotropic heat trans- 
port i n the envelope was studied in a pioneering pa- 
per by iGreenstein fc Hartkel (|1983l ) . The observational 
consequences of these models were anal yzed for th e 
pulsars Vela and Geminga among others (|Pagd fl995). 
iPotekhin fc Yakovlev! (|2001f ) calculated the angular dis- 
tribution of temperatures in magnetized envelopes taking 
into account the quantizing effect of the magnetic field on 
the electrons, and the suppression (enhancement) of the 
electron thermal conductivity in the direction perpendic- 
ular (parallel) to magnetic field lines. Nevertheless, the 
anisotropy generated in the envelope is not strong enough 
to be consistent with the observed thermal distribution 
of some isolated NSs, and it should be originated deeper 
in the NS crust. The understanding of kinetic properties 
of matter in NS crusts and envelopes has also been re- 
cently improved, with sp ecial attention received by th e 
role of ions and phonons (jChugunov and Haensell feOO?). 
which can be relevant at low temperatures and densities. 
In addition, the effect of impurities on the heat conduc- 
tion in a non-perfect lattice is also an open problem that 
must be considered in the near future. 

More recently, crustal confined magnetic fields were 
considered to be responsible for the surface thermal 
anisotropy observed in some isolated NSs. Tempera- 
ture distributions in the crust were obtained as station- 
ary sol utions of the d i ffusio n equation with axial sym- 
metry ([Geppert et al.1 l2004f) . The approach assumed 
an isothermal core and a magnetized envelope as inner 
and outer boundary condition, respectively. The results 
showed important deviations from the crust isothermal 
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case for crustal confined magnetic fields with strengths 
B > 10 13 G and temperature T < 10 s K. Same 
conclusions have been obtained considering not only 
poloid al but also toroida l components of the magneti c 
field (jPerez-Azorm et all l2006at iGeppert et all I2006T ). 
These models succeeded in explaining simultaneously the 
observed x-ray spectrum, the optical excess, the pulsed 
fraction, and other spectra l features of some isolated NS, 
such as RX J0720.4-3125 (jPerez-Azorm et al.|[2006bh . 

Although former studies about anisotropic temper- 
ature distributions on the cooling history of NSs 
(iShibanov fc YakovlevlH99aiPotekhm fc Yakovlevll2T)0lh 
provided very useful information, detailed investigations 
of heat transport in the non spherical case for different 
magnetic fie ld geometries have n ot been available until 
recently. In lAguilera et al] (|2007f ) we have revisited the 
cooling of NSs combining the insulating effect of strong 
non radial fields with the additional source of heating due 
the Ohmic dissipation of the magnetic field in the crustal 
region. We have shown that, during the neutrino cooling 
era and the early stages of the photon cooling era, the 
thermal evolution is coupled to the magnetic field evolu- 
tion, and both processes (cooling and magnetic field dif- 
fusion) proceed on a similar timescale (~ 10 6 yr). The 
energy released by magnetic field decay (Joule heating) 
in the crust is an important heat source that modifies 
or even controls the thermal evolution of a NS. There 
is indeed obse r vation al evidence of this fact. As shown 
in iPons et al.l (|2007f ). there is a strong correlation be- 
tween the inferred magnetic field and the surface tem- 
perature in a wide range of magnetic fields: from mag- 
netars (B > 10 14 G), through radio-quiet isolated NSs 
(B ~ 10 13 G) down to some ordinary pulsars (B < 10 13 
G). The main conclusion is that, rather independently 
from the stellar structure and the matter composition, 
the correlation can be explained by heating from dissi- 
pation of currents in the crust on a timescale of ~ 10 6 
yrs. 

In order to investigate if there is observational evidence 
that supports our models for the crustal field evolution, 
we have compared theoretical cooling curves including 
magnetic field decay with a sample of NSs with reported 
thermal emission. To obtain the cooling curves we have 
performed two-dimensional simulations by solving the 
energy balance equation that describes the thermal evo- 
lution of a NS 
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where C v is the specific heat per unit volume, Q v are en- 
ergy losses by ^-emission, Qj the energy gains by Joule 
heating, and k is the thermal conductivity tensor, in gen- 
eral anisotropic in presence of a magnetic field. In this 
equation we have omitted relativistic factors for simplic- 
ity. A detailed description of the formalism, the code, 
and results can be found in lAguilera et alj (|2007f ). The 
geometry of the magnetic field is fixed during the evo- 
lution. As a phenomenological description of the field 
decay, we have assumed the following law 



B = B Q 



exp(-t/r hm) 



1 + (TOhm/THall)(l - exp (-t/r hm)) 
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TABLE 1 

Properties of neutron stars with reported thermal 
emission. ordered by decreasing magnetic field strength. 
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where B is the magnetic field at the pole, Bo its ini- 
tial value, Tohm is the Ohmic characteristic time, and 



from Pulsar Online Catalogue 6 . 

"http:/ /www. physics. mcgill.ca/ pulsar/magnetar/main.html 
'http: / / www. atnf . csiro.au / research / pulsar / psreat / 



THaii the typical timescale of the fast, initial Hall stage. 
In the early evolution, when t <C Tohm, we have B ~ 
B (l + i/rHaii) -1 while for late stages, when t > rohm, 
B ~ Bq exp(— t/rohm)- This simple law reproduces 
qualitatively the result s from more complex simulations 
([Pons fc Geppertl l2007f ) and facilitates the implementa- 
tion of field decay in the cooling of NSs for different 
Ohmic and Hall timescales. 

Our sample of objects includes magnetars, isolated 
radio-quiet NSs, and radio-pulsars, and it covers about 
three orders of magnitude in magnetic field strength. Be- 
fore we discuss our results several caveats should be men- 
tioned. For most NSs, B is estimated by assuming that 
the lose of angular momentum is entirely due to dipo- 
lar radiation. If the magnetic field is constant along the 



The impact of magnetic field on NS cooling 



3 



evolution, the dipolar component can be estimated by 
B d = 3.2 x f0 19 (PP) 1 / 2 G, where P is the spin period in 

seconds, and P is its time derivative. Alternatively, for 
a few radio-quiet isolated NSs, B is estimated assuming 
observed x-ray absorption features are proton cyclotron 
lines. This association is still controversial and, if true, it 
should be read as an estimate of the surface field, which 
is usually larger than the external dipolar component. 
In order to work with a more homogeneous sample, we 
have included in our study only those objects for which 
P is available. The quoted magnetic fields are always the 
spin-down estimate for the dipolar component. 

Ages are also subject to a large uncertainty. If the 
birth spin rate far exceeds the present spin rate and B 
is considered constant, the age can be estimated by the 
spin-down age (t s d = P/2P). In the cases that an in- 
dependent estimate is available (e.g. a kinematic age), 
the spin-down age does not necessarily coincide with the 
other estimation. If one considers magnetic field de- 
cay, however, t s d seriously overestimates the true age (t) . 
Correcting the spin-down age to include the temporal 
variation of the magnetic field may help to reconcile the 
observed discrepancy between the spin-down ages and in- 
dependent measures of the ages of some of the objects we 
have included in this analysis. This makes the compar- 
ison with observations more meaningful. Other sources 
of error do probably exist but it is not the purpose of 
this letter to discuss each observation separately. 

As for the temperatures, in general there is no clear 
way to define error-bars. We have adopted the range of 
values found in the literature, and we also indicate those 
cases in which the estimate should be interpreted as an 
upper limit, more than a measure. This is the case of 
the Crab pulsar, or even of some magnetars, which show 
large variations in the flux in the soft x-ray band on a 
timescale of a few years, indicating that the real ther- 
mal component may be that measured during quiescence 
and that the x-ray luminosity during their active periods 
might be a result of magnetosphe r ic act ivity. Following 
the criteria adopted in lPage et alj (|2004l ). most reported 
temperatures are blackbody temperatures, except for low 
field radio-pulsars in which the blackbody fit results in 
an unrealistic small radius of the NS. In these four cases, 
we took the temperature consistent with Hydrogen at- 
mospheres. We include a list of the sources considered 
in Table 1, with the corresponding references. 

In Fig. [TJ we show two types of theoretical curves: 
temperature versus true age (according to the simula- 
tion) and temperature versus the spin-down age consis- 
tent with the assumptions made on each model about the 
magnetic field evolution. We have sampled the objects in 
three groups according to their measured magnetic field: 
high field NSs with Bd > 10 14 G, intermediate field NSs 
with 10 13 G< B d < 10 14 G and low field NSs for which 
Bd < 10 13 G. We found that this three samples could be 
explained qualitatively by cooling curves in three differ- 
ent regimes: high, intermediate and low magnetized NSs, 
in all cases with Joule heating included. Each of these 
regimes is represented by a set of curves with a given or- 
der of magnitude of the initial field, Bq. We have taken 
Tohm = 10 6 yr and r Ha n = 10 3 yr/B ,i5, where B ,15 is B 
in units of 10 15 G. The de pendence of the result s on the 
decay rates is discussed in lAguilera et al.l ()2007l ) . These 
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Fig. 1. — Cooling curves. Surface temperature at the pole as a 
function of t (upper panel) and t sd (lower panel) for different Bo- 
The magnetar region (red) contains curves for NSs born with Bo > 
5 x 10 14 G. Intermediate field NSs (5 x 10 13 G< B < 5 x 10 14 G) 
evolve within the brown area and the blue region corresponds to 
NSs with Bo < 5 X 10 13 G. Observations: squares for sources with 
B d > 10 14 G, triangles with 10 13 G < B d < I0 14 G, and circles with 
B d < I0 13 G. Open circles are temperatures obtained from fits to 
Hydrogen atmospheres. 



three regions are depicted in Fig. [T] as colored zones. For 
comparison, the cooling curve corresponding to a non- 
magnetized NS is represented by a dashed line. 

Focusing in the high field (red) region, we see that 
the effect of the magnetic field is visible from the very 
beginning of the NS evolution. The initial equilibrium 
temperature reached in a non-magnetized model may be 
increased up to a factor of 5 and it is kept nearly con- 
stant for a much longer time, up to I0 4 years. The effect 
of Joule heating is very significant and may help to un- 
derstand the high temperatures observed in magnetars 
(jKaspil l2007t ). although other physical processes could 
contribute as well. For instance, the initially higher tem- 
peratures result in higher electrical resistivity, therefore 
accelerating the magnetic field dissipation in an early 
epoch. We have not yet considered the consistent tem- 
perature dependence of the resistivity, in combination 
with the evolution of the magnetic field geometry. A 
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fully coupled magneto-thermal evolution code should be 
employed for this purpose (work in progress). The ob- 
served temperatures of radio-quiet, isolated NSs could be 
explained either if they are old NS (« 10 6 years) born 
as magnetars with Bo = 10 14 -10 15 G, or by middle age 
NSs born with fields in the range of Bq = 10 13 -10 14 G, as 
plotted in the intermediate brown zone of Fig. [TJ For in- 
termediate field strengths, the initial effect is not so pro- 
nounced but the star can be kept much hotter than non- 
magnetized NSs during the period from 10 to 10 6 yrs. 
Finally, the effect of the magnetic field in stars born with 
initial fields of the order of 10 13 G is moderate, but still 
visible. These would be the case of some radio-pulsars, 
for which the detection of a thermal component in their 
spectrum has been reported. For weakly magnetized NSs 
with B ~ 10 12 G, the effect of the magnetic field is 
small and they can satisfactorily be explained by non- 
magnetized models, e xcept for very old N Ss (t > 10 6 
years), as discussed in lMiralles et al.l ()1998l ). 
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Fig. 2. — Surface temperature (at the pole) as a function of the 
magnetic field. Theoretical curves are shown for different Bq and 
varying and t j^aii • The grey band corresponds to the HBL 

(Pons ct al. 2007). Observations correspond to Fig. [TJ 



The fact that the magnetic field plays a relevant role 
in the thermal evolution of neutron stars at least during 
the first million years of its life is supported but another 
observational fact. The magnetic field distribution ob- 
served in radio-pulsars is peaked in 3 x 10 12 G. Therefore, 
assuming no correlation between B and T, one would ex- 
pect to find thermal emission from neutron stars below 
and above that value in approximately the same number. 
However, almost all nearby objects for which its tem- 
perature is known have fields far exceeding that value. 
This mismatch between the pulsar field distribution and 
the thermal emitters field distribution can be easily ex- 
plained by our assumption: magnetic field decays signif- 
icantly during the first 10 6 years of a NS life and this 
effect governs its thermal evolution. As a consequence, 
NSs with larger fields have higher temperatures and cool 
slower, increasing the chances to be detected. This hy - 
pothesis has been recently outlined bv lPons et al.l (|2007l ) 
who reported a strong correlation between the surface 
temperature and the magnetic field, well approximated 



by the expression 

T s,6 - CB d,u ( 3 ) 
where T Si q is T s in units of 10 6 K, Bd,i4 is Bd in units of 
10 14 G, and C is a constant that depends on the thick- 
ness of the crust, the Ohmic dissipation timescale, and 
the ratio between the unknown internal field and the ob- 
served external dipole. For typical numbers, C ~ 10. 
This straight line in a logarithmic plot of temperature 
versus magnetic field has been named the heat balance 
line (HBL). 

To contrast our results with this hypothesis based on 
a simple energetics argument, we plot in Fig. [2] the vari- 
ation of the polar surface temperature T s (usually iden- 
tified with the hot component of the thermal spectrum) 
as a function of B for different theoretical curves, and we 
compare them with the observational data from Table 1. 
In this T-B diagram the thermal history of a NS proceeds 
as follows. A NS begins its life high on the figure with 
some initial field Bq. As it cools it moves vertically down- 
ward, until decay of its field provides an energy source 
able to counterbalance the thermal loses. This causes the 
trajectory to bend to the left. Then it continues mov- 
ing down but the temperature and the field evolution 
are coupled since heating produced by magnetic field de- 
pends on the field strength and the cooling mechanisms 
depend on the temperature. For young magnetars the 
initial stage is very short because the high field prevents 
the star to cool down further and the surface temperature 
remains relatively high. A further evolution proceeds 
only when the field is also decaying. For low magnetic 
field stars the drop in the initial temperature is more 
pronounced and they reach the equilibrium region much 
later. The region in the diagram where all cooling curves 
tend to converge i s in agreement with the HBL defined in 
iPons et ail (|2007f ) (indeed it is a band, since not all NSs 
are necessarily identical). The grey band corresponds to 
the choice C — 5-50. 

Since the initial vertical trajectory is very fast, NSs do 
not spend much time in that region of the diagram. A 
clustering of sufficiently old objects in the grey region is 
expected. Low magnetic field stars reach undetectable 
low temperatures more quickly, thus they are more dif- 
ficult to detect than highly magnetized NSs. In the late 
photon cooling era, since the main mechanism to radiate 
energy is photons from the star surface, one expects that 
the correlation given by Eq. ([3]) is more clear. Although 
it was a very simplified approximation, we show that our 
evolution curves approach the HBL slope asymptotically. 
Finally, it is worth to notice that some objects are crossed 
by curves belonging to different initial Bo, so there is no 
unique determination of the initial conditions. 

As a result of this investigation we conclude that 
consistent magneto-thermal evolution simulations are 
needed before we can disentangle the properties of the 
interiors of neutron stars by studying their cooling his- 
tory. A first step in that direction has been given in 
this work, showing that the effect of magnetic field de- 
cay in the highly resistive crust (as opposed to the highly 
conductive core) could be very large. It certainly has a 
significant impact on the thermal evolution of stars with 
B > 10 13 G. The thermal and magnetic evolution of neu- 
tron stars is (at least) a two parameter space, in which 
the evolutionary tracks of NSs born with different initial 
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conditions (the mass of the object, the initial magnetic 
field) cannot be properly described by a unique tempera- 
ture versus age curve. It becomes clear that the minimal 
model that hopefully will be used to understand neutron 
stars needs to include the magnetic field structure and 
evolution. 
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